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Abstract 

This paper presents a theory for quantum light propagation in a single-mode fiber which 
includes the effects of the Kerr nonlinearity, group-velocity dispersion, and linear loss. The 
theory reproduces the results of classical self-phase modulation, quantum four-wave mixing, 
and classical soliton physics, within their respective regions of validity. It demonstrates the 
crucial role played by the Kerr-effect material time constant, in limiting the quantum phase 
shifts caused by the broadband zero-point fluctuations that accompany any quantized input 
field. Operator moment equations— approximated, numerically, via a terminated cumulant 
expansion — are used to obtain results for homodyne-measurement noise spectra when dis- 
persion is negligible. More complicated forms of these equations can be used to incorporate 
dispersion into the noise calculations. 


1 Introduction 

Optical fibers have long been considered for the generation of squeezed-state light, starting with the 
pioneering work of Levenson and coworkers in the mid 1980’s, who observed 0.58 dB of continuous- 
wave (cw) squeezing [1], to recent measurements exhibiting over 5dB of short-pulse squeezing [2]. 
In this paper we present a theory for quantum light propagation in single-mode optical fiber. Our 
development, which includes the effects of the Kerr nonlinearity, group- velocity dispersion, and 
linear loss, is guided by two overarching principles: the theory must include all relevant prior re- 
sults, both classical and quantum mechanical, and, within reasonable limits, it must accommodate 
arbitrary input states. This theory [3] is an extension of our prior work on quantum propagation 
in a dispersionless, lossless, Kerr medium [4], In that earlier study it was shown that a material 
time constant is crucial to a correct description of quantum nonlinear phase shifts beyond the 
four-wave mixing regime, a conclusion similar to that reached by Blow and coworkers [5]. The use 
of a finite Kerr-effect time constant is retained in the current treatment of lossy, dispersive fiber. 

This paper is organized as follows. In Section 2 we review quantum propagation in a Kerr 
medium, concentrating on the structure of the theory and the necessity of a finite Kerr-effect time 
constant. We recount the principal result of [4], namely the limits on squeezed-state generation 
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in lossless, dispersionless fiber. Section 3 expands this theory to incorporate dispersion and linear 
loss via a split-step approach. Section 4 introduces the terminated cumulant expansion (TCE) 
as a technique for closing the infinite chain of coupled moment equations generated by the full 
theory. Using the TCE, the limits on quadrature-noise squeezing in lossy, dispersionless fiber 
are quantified. Finally, in Section 5, we discuss the relationship of our approach to other quan- 
tum propagation theories for single-mode fiber, focusing on the necessity of the Kerr-effect time 
constant. 


2 Quantum Self-Phase Modulation 

Our attention, in this section, is restricted to a linearly polarized field propagating in a lossless, 
dispersionless, single-mode fiber that exhibits the Kerr nonlinearity. Classically, the refractive 
index in this fiber can be written as follows, 

n(z,t) = n 0 + ^^-\E(z,t)\ 2 , (1) 

where z is the axial coordinate along the fiber, t is time, no is the linear refractive index, n 2 is the 
Kerr coefficient, A is the fiber’s effective cross-sectional core area, and E(z, t ) is the normalized 
complex envelope of the singl e-mode field within the fiber. The field normalization we employ is 
such that E(z , t) has units yphotons/sec. Note that the introduction of the photon energy, hu, 
is strictly a convenience at this classical stage; it has no quantum significance as yet. 
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Fig. 1. Schematic configuration for Kerr-effect 
propagation in lossless, dispersionless fiber. 


The classical propagation problem in lossless, dispersionless, Kerr-effect fiber is sketched in 
Fig. 1. The fiber is excited, at z = 0, by an input field E IN (t) that launches a -l-z-going wave 
E(z,t ) satisfying E(Q, t) = Eis(t). In a reference frame moving at the group velocity, v g , the 
complex field envelope within the fiber satisfies the following differential equation [6], 


dE(z, t') 
d~z 


inE*{z , t')E(z, t')E(z , t'), for z > 0, 


where t' = t — z/v g is the retarded time, 


k = 


27rn 2 hu; 
AX ’ 


(2) 

(3) 
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is the nonlineax phase shift per unit length per unit photon flux, and A = 27 tc/o; is the center 
wavelength of the light. The intensity is a constant of motion for Eq. 2, so it is easily shown that 

E(z, t) = exp[i/c.z£*(0, t')E(0, t')\ E( 0, t'), for z > 0. (4) 

Using the initial excitation condition that specifies E{ 0, t 1 ), and directing our attention to Eour(t ) = 
E(L, t), the field coupled out of the fiber at z = L, we obtain the classical input-output relation 
for a length L lossless, dispersionless fiber exhibiting the Kerr nonlinearity, namely, 

Eour{t) = exp[irE} N (t)E IN (t)] E IN (t), (5) 

where 

r = kL , (6) 

is the nonlinear phase shift per unit photon flux and, for convenience, we have dropped the L/v g 
group delay. 

Two well known results follow directly from Eq. 5: spectral broadening through self-phase 
modulation (SPM), and four- wave mixing (FWM). As an immediate consequence of Eq. 5, we 
see that an optical pulse propagating through the fiber acquires a time- varying phase shift, which 
is proportional to the pulse’s intensity. The derivative of this time-varying phase constitutes an 
intensity-dependent instantaneous-frequency variation, implying, for sufficiently intense pulses or 
long fibers, significant spectral broadening. On the other hand, when the input field comprises a 
strong monochromatic pump, Ef N , at frequency u plus weak sidebands at frequencies o;±fl, FWM 
couples the sidebands. The input-output relation for classical FWM can be obtained, from Eq. 5, 
by replacing the exponential term with its two-term Taylor series approximation, and assuming 
that all nonlinear phase shifts other than the pump x pump term are small: 

Eour(t ) = exp(fr|Ef w | 2 ) [(1 + ir\Ef N \ 2 )E IN (t) + ir(Ef N ) 2 E* IN (t)] , (7) 

where E(t) = E{t) — E p for the input and output fields. 

In the quantum theory of the Kerr interaction the photon-units complex field envelope, E(z , t), 
becomes a photon-units field operator, E(z,t). The input field operator — Emit) = E(0,t) — is a 
single-spatial mode, multitemporal-mode free field, and hence must satisfy the following 5-function 
commutator rule [7], 

[£/*(*), £/ W (u)] = S(t - u ), (8) 

for photon-units field operators. The output field operator — Eour(t ) = E(L, t) — is also a single- 
mode free field, whose commutator must therefore mimic Eq. 8. This requirement is automatically 
met by quantizing Eq. 2 as follows, 

— — - = ikE\z, t')E{z, t')E(z, t'), for z > 0. (9) 

Quantum FWM emerges from Eq. 9, in a manner that ensures commutator preservation, by 
decomposing the input field operator, E IN {t), into a c-number pump, Ef N , plus a 5-function 

commutator field operator, Ej^{t). Then, assuming that all nonlineax phase shifts other than the 
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pump x pump term are small, Eq. 9 can be linearized and solved to yield the quantum version of 
Eq. 7, namely, 

Eour{t) = e x p(tr|is^| 2 ) (1 + ir\E? N \ 2 )EiN(t) + ir(Ef N ) 2 E IN (t) . (10) 

Equation 10 is commutator preserving, does not require a Kerr-effect time constant, is consis- 
tent with the classical FWM theory, and agrees with experimental results in fiber squeezed-state 
generation [1]. 

We find that it is not possible to treat quantum SPM from Eq. 9, in a way that recovers the 
classical limit, without a modification of the theory. To highlight the failing of instantaneous- 
interaction quantum SPM, let us calculate the mean output-field for a coherent-state input field 


[4], By direct substitution, 

E(z,t) = exp [iKZ&{0, t')E(0, t')] E( 0,f'), for z > 0, (11) 

can be shown to be the solution to Eq. 9, leading to 

Eour{t ) = exp ^irE\ N (t)Eipf(t)^ Ei^(t), (12) 

as the quantum version of Eq. 5 under an instantaneous Kerr-interaction model. The output mean 
field, 

Eour{t ) = (EiN(t)\EouT(t)\EiN{t)) 

= (E IN (t)\exp^irE] N (t)E IN {t)^E IN (t)\E IN (t)). (13) 

for a coherent-state input field reduces to 

EouAt) = {Emit ) | exp [irE\ N (t)E IN (t)} \E IN (t))E IN (t), (14) 

because a coherent-state input field obeys the eigenfunction relation 

E/iv(t)\EiN(t)) — Ei^(t)\EiN(t)), for all t. (15) 

For the rest of the calculation we employ a limiting argument, 

EouAt) = hrn ^(i?/jv(f)| exp (ir/T) ^dr ■E|j V (t).E/jv(t)J | Ej^^t))^ Eji^(t) (16) 

= lim ^exp (exp(zr/T) - 1) £ ^dr \E in (t)\ 2 J^ E IN (t) (17) 

= E IN (t), (18) 


where we have exploited the characteristic function for coherent-state photon counting [8] to obtain 
Eq. 17. 

Equation 18 predicts that the instantaneous Kerr nonlinearity has no effect whatsoever on the 
mean field of a coherent-state input, regardless of how large the classical, peak nonlinear phase shift 
becomes! This result contradicts the classical theory, Eq. 5, and, more importantly, it contradicts 
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experiments which have confirmed the classical spectral-broadening predictions. In retrospect, 
the failure of Eq. 9 to reproduce SPM in the appropriate limit should not be too surprising as it 
arises from applying an instantaneous, spectral-broadening nonlinearity to a quantized field whose 
vacuum-state fluctuations extend to infinite bandwidth. 

The classical SPM result can be recovered from our quantum mean-field calculation by intro- 
ducing a phenomenological Kerr-effect time constant, tk, as a lower-bound for T in Eq. 16. In 
particular, with T = r K > 0 and rjr K < 1, Eq. 17 becomes 


Eour{t ) ~ exp 


(ir/r K ) [ < 

L h-T K 


(1t\E in (t){ 


E IN {t) 


exp[ir|.E'/7v(t)| 2 ] E IN (t), 


(19) 


whenever tk is smaller than any classical time scale present in the field. 

To incorporate t k into our quantum propagation theory for the lossless, dispersionless, Kerr 
effect in single-mode fiber, we employ the following partial mode expansion for the input-field 
operator: 

OC 

Ew(t) = a ! n N t(t - nr K ), (20) 

71= — OO 

where { a^ N oo<n<oo}isa set annihilation operators, and 


m = 


i 

V*' 


for 0 < t < tk, 


{ 0, otherwise, 


( 21 ) 


is a c-number, T/f -second duration, rectangular pulse. Suppressing the L/u 9 -second group delay, 
the same expansion applies to the output field operator, E OUT {t ), with the input annihilation 
operators replaced by the output annihilation operators, { a° UT : — oo < n < oo }. These output 
annihilation operators are related to their corresponding inputs by 


&n UT = exp [(ir/rjrja^o^ a™, for -oo < n < oo, 


( 22 ) 


which should be compared with the instantaneous quantum Kerr model, Eq. 12. 

Equations 20-22 comprise a coarse-grained time model for quantum Kerr-effect propagation 
in lossless, dispersionless, single-mode fiber. Introducing the Kerr time-constant as a simple phe- 
nomenological parameter, tk, is arguably a distasteful ad hoc procedure, but the known existence 
of finite Kerr response times plus the inconsistency of an instantaneous quantum Kerr model mit- 
igate against this deficiency. A more important consideration with respect to Eqs. 20-22 is that 
they are fundamentally incomplete; {£(£ — r itk) ■ — oo<n<oo}isan orthonormal temporal 
mode set, but not a complete orthonormal mode set. Hence, in the coarse-grained time model, 
the input field commutator obeys 


[E IN (t),E\ N (u)] = 6 TK (t,u) 

{ ■f£, if ( n — 1 )tk <t,u< utk , for n an integer, 

0, otherwise. 


(23) 
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This expression is not the 6-function commutator that a photon-units field operator should possess. 
However, as will be made clear below, it can be a reasonable approximation thereto. Moreover, 
from Eqs. 22 and 23, we see that the coarse-grained time quantum theory is consistent in that it 
preserves 6 Tk -function commutators — we have that 

Eour(t), Eq UT (u)} = 6 TK (t,u). (24) 

Given the approximate nature of the coarse-grained model’s field commutator [4], what con- 
ditions suffice to ensure the accuracy of this model’s moment predictions? Clearly, this model is 
a poor approximation for field dynamics occurring on time scales comparable to t K) or, equiv- 
alently, over band widths of the order of r^ 1 . We expect t k to fail in the interval 1-100 fs— to 
match the known time scales of n 2 — so our quantum input field, E IN (t ), must not have excited 
(non-vacuum state) modes at these time or bandwidth scales. Furthermore, because SPM causes 
spectral broadening, we must check that this same time/bandwidth scale condition is satisfied 
by the output field operator, Eourit )• Finally, we must limit any photodetection measurements 
that we make on the output field to be insensitive to i?oi/T (^)~behavior at tk time scales or 1 
bandwidths. When all three of these conditions apply, we believe the coarse-grained time quan- 
tum Kerr model — with a proper value for t k — should capture the full behavior of quantum Kerr 
propagation in lossless, dispersionless, single- mode fiber. 

The output mean-field derivation — with a Kerr time-constant t k — is virtually identical to the 
analysis of the instantaneous-interaction case. We find that 

EouT^t) = ( Eiw(t)\EouT{t)\EiN{t )) = exp^ii?|F//^(f)| j Ei^(t), (25) 

where we have defined the complex number, R, according to 

iR = (exp(# 9 ) - 1 ]t k , for 4> q = rjr K . (26) 

Equation 25 bears a striking similarity to the classical formula, Eq. 5. Indeed, the two would be 
identical were R — r to prevail. 

We know that r is the classical nonlinear phase shift per unit photon flux. Because t k is the 
time duration of a single, rectangular- pulse mode in our coarse-grained time model, we have that 
(j)q is numerically the classical phase shift produced by one photon in such a mode. Equation 25 
applies to a coherent-state input field— which has vacuum-fluctuation noise in all its modes— so 
it is fair to regard (p q , physically , as the quantum nonlinear phase shift that is due to the vacuum 
fluctuations of a single mode. 

The predictions of the quantum and classical theories are nearly coincident when we have 
r « R e{R) and |Im(R)| < 1. These conditions hold for <p q «: 1, as can be seen from Fig. 2, 
where we have plotted R e(R)/r and Im(R)/r vs. (j) q . Assuming t k = lfs, a 1km fused silica 
fiber will have 4> q « 10“ 3 rad; from Fig. 2 we conclude that Eq. 25, the quantum coherent-state 
mean field, will then be in excellent agreement with Eq. 5, the classical mean field result. For 
<p q > 10" 1 rad, there is a pronounced divergence between R e(R) and r, caused by the intrinsic 
periodicity of R. Viewed as a function of <p q , Eq. 26 shows that R is periodic with period 2n. 
Physically, this periodicity constitutes a quantum state-recurrence for Kerr-effect propagation in 
lossless, dispersionless fiber and occurs when L = L q = 2ittk/k', for tk = 1 fs and fused-silica 
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fiber this implies L q ~ 10 4 km, an experimentally inaccessible value for our assumptions of zero 
loss and zero dispersion. 


i 



A 

2k 


Fig. 2. Logarithmic plots of Re(R ) / rand Im(R ) / r 
vs. q. The coherent-state mean field of coarse- 
grained time SPM reduces to the classical SPM 
result when R « r, q is the quantum nonlinear 
pfiase shift, i.e., the Kerr-induced phase shift of 
one photon per K -second mode. 
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Fig. 3. Minimum and maximum low-frequency, 
homodyne-measurement noise spectra for a cw, 
coherent-state input vs. fiber length. The solid 
curves are the coarse-grained time SPM the- 
ory; the dashed curves are the instantaneous- 
interaction FWM theory. The parameter values 
employed are: = 1.06 pm; A = 3 56 x 10' 11 m 2 ; 

n 2 = 32 x 10“ 20 m 2 /W; P IN = 1W; and K = 
lfe. 


Expressions for the quantum output moments, up to the second order, have been derived for 
Gaussian-state inputs [4], We note that these moments yield the correct results in the appro- 
priate limits. Specifically, when <fi q 1 they are in agreement with the quantum FWM results. 
Furthermore, if we let R — ► r and S Tk (t, u) — > 0 and use classical covariances in lieu of quan- 
tum covariances throughout, e.g., (E] N (t)Ei^(u)) — > Ej N (t)E[^(u), etc., the resulting equations 
agree with classical stochastic SPM results. We note that none of these classical formulas depend 
on the Kerr-effect time constant, Tr. This is fully consistent with our assumption that the spec- 
trum of the classical input excitation and the implied output-field spectrum are both narrower 
than t^ 1 . Under these circumstances we expect — in a classical theory — that the Kerr interaction 
is effectively instantaneous. 

The implications of the full theory are readily ascertained by examining the output field’s 
homodyne-detection statistics. An ideal homodyne measurement on Eour(t) yields a photocurrent 
whose statistics are proportional to those of the following abstract quantum measurement: 

£^(f) = Eour{t )e*^ + E,Q UT {t)e ~ **, (27) 

where 4> is the local-oscillator (LO) phase. As is conventionally done in cw squeezing experiments, 
we shall focus on the minimum and maximum values of the homodyne-noise spectra as the LO 
phase is varied. 
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For a coherent-state input, both the instantaneous-interaction FWM and the coarse-grained 
time SPM theories imply frequency-independent homodyne spectra out to frequencies comparable 
to r^ 1 . In Fig. 3 we have plotted 5' mm (0) and ^^(O), vs. fiber length L, for both theories. These 
curves assume n<i = 3.2 x 10 _ 2 O m 2 /W, A — 3.56 x 10 -11 m 2 , A — 1.06 //m, Pjn = tuj\Ei^\ 2 = 
1 W, and, for the coarse-grained time case, t k — 1 fs. Quantum FWM predicts complete noise 
suppression for one quadrature component, 5 m j n ( 0 ) — * 0 , as the fiber length becomes infinite. 
Moreover, this transpires at minimum uncertainty product, -S' m i n (0)5' mai (0) = 1 . These FWM 
features are clearly evident is Fig. 3; for long enough fibers, however, they are at odds with our 
coarse-grained time SPM theory. According to Fig. 3, the coarse-grained time SPM theory has 
an 5 m i n (0) which reaches a nonzero minimum at L ~ 500 m, and increases, for longer fibers, 
substantially above the coherent-state value of unity. 

The precise location of the minimum in the low noise quadrature depends on tk with higher 
values of tk extending the region of agreement between the instantaneous interaction FWM and 
the coarse-grained time SPM theories. For various technical reasons, e.g., guided acoustic- wave 
Brillouin scatter, it is not feasible to perform such an experiment using a cw field. Nevertheless, 
the physical conditions corresponding to the minimum in Fig. 3, namely 1 Watt input power and 
roughly a 500-meter-long fiber, seem well within the realm of possibility for pulsed experiments. 
In fact, modeling Shelby’s soliton-squeezing experiments [9] by using his peak intensity as the 
intensity in our cw theory, we find that conditions are right for a r^-dependent deviation from 
quantum FWM. However, accurate analysis of a short-pulse experiment — especially one based 
on solitons — must surely account for group- velocity dispersion. In addition, a realistic quantum 
propagation theory for long fibers should address linear loss. 


3 Inclusion of Dispersion and Loss 

A split-step configuration for incorporating group-velocity dispersion (GVD) and linear loss (LL) 
into quantum propagation analysis for a Kerr-nonlinear fiber is shown schematically in Fig. 4. An 
infinitesimal length of fiber at 2 is divided into two sub-segments. The first sub-segment exhibits 
only the Kerr nonlinearity and linear loss; the second has neither Kerr effect nor loss, but suffers 
from dispersion. 

Let E(z,t ) be the coarse-grained time, + 2 -going, photon-units field operator within the fiber, 
i.e., as before we have the mode expansion 

OO 

E(z,t)= a n (z)£(t-nT K ). (28) 

71 — — 00 

Over an infinitesimal, 52-meter-long fiber segment, the split-step procedure leads to the 2 -to- 
(2 + 8z) annihilation operator transformation, 

OO 

a n (z + 8z) — ^2 M n — m ; ^ 2 ]a m ( 2 ), for —00 < n < 00 , (29) 

771 = — OO 


where 

a n (z ) = exp^(/c52/rK)aJ,(2)a n (-2 : ) + ( 0 / 2 ) 62 ] a n (z ) + \fabzb n {z ), for —00 < n < 00 , (30) 
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characterizes SPM and LL alone over the dispersionless sub-segment, and 

— for n = — 1,1, 

2 rfc 

h[n-8z] = \ i + ifysz, for n = 0, ( 31 ) 

t k 

L 0, otherwise, 

is a discrete-time impulse response accounting for the dispersive sub-segment. In Eqs. 28-31, we 
have suppressed the group- velocity delay, and we have introduced the fiber’s dispersion coefficient, 
/? 2 , and its power-attenuation coefficient, a. More importantly, with the inclusion of a Langevin 
noise-operator b n (z ) required by the presence of the LL these equations preserve the b TK (£, u)- 
function commutator of a coarse-grained time input field operator, Eiff(t). So, if we drive the fiber 
at z = 0 with such an input — forcing E( 0, t) = E IN (t) — and then iteratively apply Eqs. 28-31, we 
will arrive at an output field operator, Eout{ 1) = E(L,t ) at z = L with a proper coarse-grained 
time commutator. Furthermore, because the unitary operators that transform the (d n (z)} into 
the (a n (z)} are known, as are the operators for changing the (a n (z)} into the {a n (z + Sz)}, we 
can — in principle — calculate all the measurement statistics for Eout{1), given any state of E[^(t). 


Input Plane Output Plane 

<z=0) z (z=L) 



Fig. 4. Schematic split-step configuration for in- 
clusion of group- velocity dispersion and linear loss 
into the coarse-grained time SPM theory. 


We note that the only aspect missed in this split-step approach is a term proportional to the 
commutator of the Hamiltonians governing each step. It can be shown that this commutator is 
finite and its contribution goes to zero in Sz 2 ; hence, in the limit of Sz — > 0 this is an exact 
theory. Taking the limit of Sz — * 0, the differential equation for the mode operator is 

<9 3 

— d n (z) = -(a/2 )an(z) + i{K/T K )al{z)a 2 n (z) + ^-[2 a n (z) - a n _ x (z) - a n+l (z)] 

+ \/a£> n (z), for 0 < z < L, — oo < n < oo. (32) 

The coarse-grained time cannot be suppressed — it is essential to preventing the mean-field contra- 
diction exhibited in Sect. 3.1 — but it can be hidden. Returning to field operator notation, Eq. 32 
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can be recast as 

§- z E(z,t) = iK&(z,t)E(z^)E(z,t)-ify-^E(z,t)-%E(z,t) + J5F(z,t) (33) 

with the obvious, implied definitions for the -approximation to the second partial derivative 
with respect to time and for the Langevin noise field-operator. 

Equation 33 is extraordinarily appealing. Converting it, naively, to continuous-time classical 
form by merely dispensing with the operator carets, dropping the Langevin noise source, and the 
tk subscripts, we obtain 

^E(z, t) = »«£•(*, t)E(z,t)E(z,t)-i^E(z,t)-?E(z,t), (34) 

the well-known starting point for the classical theory of solitons in fiber [6]. 

4 Moment Propagation and the Terminated Cumulant 
Expansion 

Although it is possible to calculate the exact state transformation for the preceding quantum 
theory, it is a daunting numerical task in almost all cases of interest. Thus, we have elected to 
follow a much simpler and restrictive course — moment propagation. Taking the expectation value 
of Eq. 33 we find that the mean field develops according to 

£(E(z,t)) = m(E'(z,t)E(z,t)E(z,t)) - i^^(E(z,t)) - §<£(*, f)>, (35) 

which illustrates the fundamental problem of moment propagation — the Kerr nonlinearity couples 
each moment’s differential equation to those of higher order. For example, in the single time case, 
the differential equation for the moment (E^ k (t, z)E l (t, z )) includes terms containing the moment 
(£t k ^(t,z)E l+1 {t, 2 )), leading to an infinite progression of coupled differential equations. 

This infinite linkage of moment equations can be broken, in an approximate way, through 
a terminated cumulant expansion (TCE). (A brief review of cumulants — for classical random 
variables — is found in the Appendix. Here, we rely on normally-ordered quantum cumulants.) In 
the TCE- A expansion, all normally-ordered quantum-field cumulants beyond the A-th order are 
set to zero: 

«n[#(*.*)r n[^( z > *)]"*)) = °> when + n k) > k - ( 3e ) 

k k k 

The TCE- A assumption provides low-order moment expressions for all field-operator moments 
beyond A-th order. For example, the third-order cumulant relation, 

((E'(z,t)E 2 (z,t))) = (E\z,t)E 2 (z,t))-2(E\z,t)E(z,t))(E(z,t)) 

- (E'(z,t))[E\z,t))+2(E'(z,t)){E(z,t))\ (37) 

affords the following explicit expression for the TCE-2 expansion: 

(Ef(z,t)E 2 (z,t)) = (EHz,t))(E 2 (z,t))+2(P(z,t)E(zJ))(^J))-2(EHzJ))(EU^t)) 2 - (38) 
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Substitution of this expression into the mean-field equation, Eq. 35, eliminates the third-order 
moment and leads, ultimately, to a closed system of differential equations for 

{(E(z, t)>, <£t ( z , t )E{z, u)), (E(z, t)E{z , u))}. (39) 

The accuracy of the TCE approximation depends upon both the initial state and its subsequent 
propagation. For a Gaussian state, such as a coherent state, all cumulants of order three or higher 
vanish. Moreover, a Gaussian state remains Gaussian under linear propagation — even if it is lossy 
and/or dispersive — so, for example, TCE-2 is exact for Gaussian-state inputs in the four-wave 
mixing limit. Higher-order TCE approximations track deviations from a Gaussian state, hence 
they should prove useful for Gaussian-state inputs even beyond the four-wave mixing regime. T his 
is demonstrated, quantitatively, in Fig. 5, where we compare coherent-state input, homodyne- 
noise output spectra for a lossless, dispersionless Kerr-effect fiber computed via TCE -K, with 
the exact solution presented in the previous section. The coarse-grained time SPM curve (exact 
solution) represents a state that is very nearly Gaussian state up to the point of its minimum-noise 
curve departs from the instantaneous-interaction FWM curve. (As noted above, coherent-state 
FWM is always a Gaussian-state case.) The TCE-2 approximation misses the mark, as it always 
represents a minimum-uncertainty Gaussian state, which is plainly a bad approximation to the 
exact solution. However, the TCE-3 approximation captures the essential nature of the exact 
solution, viz. it correctly predicts the minimum noise level and its subsequent rise to the shot- 
noise-level. As expected, the TCE-4 and TCE-5 approximations show slightly better performance, 
but the meager improvement they provide hardly justifies their added computational burden. 



Fig. 5. Minimum and maximum low-frequency, 
homodyne-measurement noise spectra vs. q /2 
( q = 2 for a fiber of length 3 25 x 10 5 km, 
k — lfs) at 1.55 pm in a dispersionless, lossless 
fiber. SPM denotes the exact calculation, FWM 
the four- wave mixing approximation, and TCE 
the terminated cumulant expansion (K = 2,3,4,5 
shown). The input field is a 1-Watt, coherent 
state; 0 dB is the coherent-state noise level. 



Fig. 6. Minimum and maximum low-frequency, 
homodyne-measurement noise spectra vs. q /2 
( q = 2 for a fiber of length 325 x 10 5 km, 
k = lfs) at 1.55 pm in a dispersionless fiber. SPM 
denotes the exact, lossless calculation, and TCE 
the terminated cumulant expansion (K =3,4,5 
shown) with 0.2dB/km linear loss. The input field 
is a 1-Watt, coherent state; OdB is the coherent- 
state noise level. 


Armed with the TCE approximation, we can obtain homodyne-noise spectra for situations in 
which the exact calculations are thwarted by moment linkage. Consider propagation in a lossy, 
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dispersionless Kerr-effect fiber. In Fig. 6 we have plotted the exact SPM result for the lossless 
case (solid curve) and the TCE solutions (dashed curves) for a fiber with a power loss coefficient 
of 0.2dB/km. The regions of overlap follow the same trends seen in Fig. 4, and predict a 4dB 
increase in the noise level over the lossless case. 

Inclusion of dispersion couples time slots and greatly increases the complexity of TCE moment 
calculations, even with a coherent-state input. If we address N time samples of the field, there 
are (19iV + 15iV 2 + 2iV 3 ) /3 complex moments, or (35 N + 30iV 2 + 4iV 3 )/3 real quantities, to track. 
Each real quantity obeys one of 23 types of differential equation, which contain anywhere from 
5 to 87 terms. For N — 100, there are thus 717,300 moments, or 1,434,500 real quantities, to 
be computed. We are working on these calculations at present, and expect to be reporting our 
results in the near future. 


5 Conclusions 

We have presented a general theory for quantum propagation of an optical field in a lossy, dispersive 
Kerr-effect fiber. Our approach leads to equations that are continuous in space, but discrete in 
time. The time granularity is set by a phenomenological Kerr-effect time constant needed to 
properly recover the known results of classical self-phase modulation. Other theories have been 
developed that describe propagation in such a fiber [11], [16], but ours is the first for which a 
material time constant has been specifically employed to temper the instantaneous interaction. 
It has been argued that the presence of dispersion provides a much more constrictive bandwidth 
limitation than t^ 1 , thereby eliminating the need for this Kerr-effect time constant [17]. We 
disagree, but in the interests of brevity, we shall confine our remarks to a few brief points. First, 
it has been shown that there is a four-wave mixing region in which dispersion enhances squeezing 
[18]; here we may expect that dispersion exacerbates the need for a finite tk to correctly determine 
the validity limit of FWM. On the other hand, if there are propagation regimes — such as soliton 
propagation — wherein dispersion renders a finite Kerr-effect time constant unnecessary, then that 
impotence should appear in our calculations, i.e., our noise results should be insensitive to the value 
we assign to r K . Note that, even with loss and dispersion, the value of t k is irrelevant to linearized 
noise analysis, and this includes the linearized noise theory of quantum solitons. Finally, the 
theory we have presented handles the case of an arbitrary field, in either the normal or anomalous 
dispersion regimes, and is more encompassing than those restricted to soliton propagation. 


Appendix 


For a real- valued classical random vector X = (Xi, X 2 , ■ ■ ■ , X n ) whose joint characteristic function 
is (f>(s) = (expos' • X)), the cumulants are defined by: 


«n*r*» 



(40) 


where the cumulant generating function is $( 5 ) = ln(0(s)). Higher order cumulants contain 
information of decreasing significance [10]. 
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